354 



Discontinuous Fluid Motion past a Bent Plane, with Special 

Reference to Aeroplane Problems. 

By G. H. Bryan, Sc.D., E'.E.S., and Egbert Jones, M.A. 

(Received November 26, 1914.) 

The following investigation is based on the Kirchhoff-Helmholtz theory of 
discontinuous stream-line motion of a fluid moving through an orifice or 
past a barrier. The problem is two-dimensional and the fluid is assumed to 
be frictionless and incompressible. 

The question has been discussed fully by Sir G-. Greenhill, F.R.S., in his 
report to the Advisory Committee for Aeronautics, and in his book on the 
Dynamics of Mechanical Flight. He considers at length the fluid motion 
past a barrier, but in what follows a bent barrier will be considered. 

G-reenhill refers to a problem similar to this, where the stream is supposed 
to divide at the bend. The present case can be reduced to that one, by 
making two parametric angles equal (P = B), 

Integrals are obtained giving the lengths of and thrusts on the planes, and 
by means of approximate methods these will be evaluated. It will also be 
shown how these integrals can be evaluated exactly, but the expressions to 
which they have been reduced are so long that they cannot be of great value 
in investigating special cases. 

The results arrived at will be used in finding the lift and drift on the surface, 
and the relative advantages of planes bent at different angles and along 
different lines (parallel to the edges) will be considered. 

Supposing the surface to be the main plane of an aeroplane, and the fluid 
to be air, the numerical results cannot be expected to be strictly quantita- 
tively correct, inasmuch as the motion is supposed two-dimensional and the 
fluid here considered is assumed to be frictionless and incompressible, which is 
by no means true for actual fluids. Consequently, instead of a region of dead 
water, separated by a surface of discontinuity from the moving ]3art, eddies 
are set up in the wake of a surface moving through air. 

Qualitative conclusions can, however, be drawn, and the effect of camber 
compared for various planes, and in this respect our results appear to agree 
with experimental results. 

§ 1. List of the Ftmdions employed. 

Although the theory of discontinuous motion is well known and explained, 
not always very lucidly, in text-books, it will be convenient to briefly 
summarise the functions and transformations employed, as follows : — 
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z = X + iy] X and y are the co-ordinates of a point in the z diagram, being 
a representation of the cross section of the planes and the fluid. 

w ^ (jy -^ iyfr] (f) velocity potential, ^fr stream function. In the w diagram 
lines parallel to the real axis represent the stream-lines. 

f = Q (dz/dw), where Q is the velocity at infinity, 
mod f = Q/q^ where q is the resultant velocity at any point, 
arg ^ = 0, the directional angle of this velocity, 
n = log^=log{Q/q)-^id, 

In the n diagram fixed straight boundaries are represented by the lines 
6 = constant parallel to the real axis, and free boundaries of surfaces of discon- 
tinuity by straight lines parallel to the imaginary axis (since by the pressure 
equation q is constant) usually by the imaginary axis itself log (Q/q) = 
corresponding to ^ = Q. 

In the case of a stream of infinite breadth impinging on a single barrier, 
and dividing, at a point P in the barrier, into parts flowing past the two sides, 
the dividing stream line in the w diagram invariably takes the form shown in 
fig. 3 (§ 2), and it is necessary to apply the transformation t = ^w in order 
to reduce the diagram to the form of fig. 4 (§ 2), the angle between PD and 
PE thus being changed from 360° to 180°. 

This transformation is, of course, a case of the Schwarz- Chris toffel trans- 
formation 

dn/dt = const (t- ay ^~^(t-hy^^~^(t---c)y^''-K..,^ 

by means of which points on the sides AB, BC, ... , A of a closed rectilinear 
polygon in the ft diagram are represented by points on the real axis in the t 
diagram ; a, h, c, etc., corresponding to the corners A, B, C, etc., of the polygon 
at which the internal angles are u, /3, 7, etc. Points within the polygon are 
represented by points in the positive half plane in t. 

As t increases through the value a say the corresponding factor of the 
modukis changes hom.{a — tYl'"~'^ to (t—a)"-!''''^, i.e. a factor ( — 1)"/'"-^ is added 
to dn/dt, so that the argument in the ft diagram changes by tt— a. 

The formula can of course be written 

dn/dt = const, {t - ayf^ (t - hff^ . . . , 

a', /3\ etc., being the external angles corresponding to a, yS, etc. 

If we assume then that the relations connecting w and t, and ft and t, are 
known the problem is solved. 

Let dia/dt = F (0, dn/dt = G (0, 

ft = log f = log (dz/dw), when Q, = 1. 

^ Forsyth, 'Theory of Functions,' § 267 (1900). 

2 F 2 
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But 



a 



G (0 dt. 



Therefore ^ = ef^ (*> ^\ i.e, dz/clw = /» W ^* . 

Therefore dz/dt = e-^^ ('> ^' ¥ (t), 

which gives a relation between the length of the plane and t 

dt, taken over the front and back of the 



dz 



For the thrust we have \p 

plane, so that if the density of the fluid be p, the thrust is 

dz 



P 
2 



il-f) 



dt 



dt, when Q = 1. 



§ 2. The Theory of Discontinuous Motion Applied to a Bent Plane. 

After having given a list of the functions employed, we proceed to use 
them to find the thrust on a bent lamina. 

Let ABC be a section of the lamina, the relative motion of the fluid being 
denoted by the arrows (fig. 1). Let a and yS be the angles AB and BC respec- 
tively make with the direction of the fluid at infinity. Let P be the point 
at which the stream divides. 




-<^^ 



. i| ^ .... w ^^^ 



Fig. 1. 

The velocity vanishes at P. The velocity also vanishes at the bend B, 
and must therefore have a maximum at some point E in BP. 
Let Q = 1 (until | 6), then 

O = log(l/q)-i'id. 

If ^ = a along PB, = ^ along BC and a— tt along PA, 
At D and E ^ = 1 and ^ = 0, and g' = at P and B. 
Consequently the XI diagram is as shown in fig. 2. 



* Lamb, ^ Hydrodynamics' § 63 (1906). The angle ABC (above, fig. 1) <7r. 
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w = cjb+i'^ is represented in fig. 3, and we assume the integration constants 
so chosen that t^ = at P. 



P 



I Jo 







T) 



xx 



E 



Fi0. 3. 



Finally, with t—p = ^w, we obtain the f diagram, as given in fig. 4, 



E at CO A 






R B 



C 



D at 00 



Fig. 4. 



The Schwarz-Ohristoffel formula is now applied to transform the O diagram 
to the t diagram. In this case it is 



do, 



const, (^—r) 



dt {t-l)){t-p)\/{t-a){t'-'ey 

where t = r, b, p, a, and c are the values of t at E, B, P, A, and C 
respectively. 

Splitting up by partial fractions and putting ^ = l/(^^— &) in one case, and 
= l/(t—p) m the other, and integrating we have 



n = Aiiog^ 



+ A^ loff \/(i-(^)(o-p) + \/(t-e)(a' 



-p) 



\ / 



Referring now to the O and t diagrams, figs. 2 and 4, it is evident that as t 
increases through P, Xi increases through tt, and as t increases through B, 
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H increases through /3— a. Now, calling y8— a = 7r/2v and neglecting the 
constant factor which occurs throughout, we may write 

n - o 1 n^ \/(^ - ^) (^ " ^) + \^i^ - g) (g - &) 

Let us now choose the scale of representation so that c—a = 1. Then for 
c'^t^ a there is a real angle U between and |-7r, such that G—t = cos^ U 
and t — a^ sin^U.* 

This gives us 

i. ^=,(t -pf = (sin^ U - sin^ P)^ = sin^ (U + P) sin^ (U - P). (3) 

ii. ^y(t'— a) {G—'b)-{-^(t'—G) (a—h) =. sin (U ^B), and so on. 

Consequently, 

n = wSin(U + P)^i.i sin(U-fB) 

^sin(U-P) 2z. ^sin(U-B) ^ ^ 

The Angles of Attack, — By the term the " angles of attack " we understand 
the angles the two parts of the bent plane make with the direction of the 
stream at infinity. 

Any particular plane has the following characteristics : P, B, and v, B 
depends on the ratio of AB to BO (fig. 1) and v on the angle between AB 
and BC. We now proceed to find an expression for P by considering the 
angles of attack. 

AB will be henceforth called the front plane, and BC the rear plane. 

Now 

a = Axlog^^(^-^^)(^7^>+/(^-:^>(^^"^> 

+ A. log >/(^"^)(^,7^)V (^--^>^^-^^> = log 1 + ^0. 

First of all, when c > ^^ > 6 it is evident that, with c > & >j^ > a, the right- 
hand side of the above equation is real, therefore 

6> = 0. 
This is along the rear plane. 

Again, when h^t'^'p, the part whose coefficient is A2 is real, whereas the 
other is imaginary and = Ailog(l/^). 

Therefore = Ailog(l/^) = -Ai7r/2. 

This is along PB. Finally, p^t^a gives us 

e = {Ai + A2)log{l/i) = ^(Ai + A2)7r/2 
along AP. 

* Cf. Greenhill ' The Dynamics of Mechanical Flight.' 
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From fig. 5 we see that if = along BC, the values of 6 along BP and 
AP are —7r/2v and — (7r + 7r/2z^) respectively, giving 

Ai = 1/v and A2 = 2, 
as before. 




Fig. 5. 



Again, for t = 00 , 



= Ai log (cos B -f^ sin B) + A2 log (cob P + '^ sin P) 

Hence = 2P4-B/z;. (5) 

This, then, is the angle of attack on the rear plane. 

We will confine ourselves to positive angles of attack. The angle of attack 

on the front plane is 

2P + B/z/-7r/2i/. (6) 

This is by assumption positive. 
Hence 2P > 7r/2^-B/^. (7) 

Hence the stream must divide beyond a certain point defined by 

2P = Itt-B/z/, 

which is the limiting position in the case of the front plane being placed 
parallel to the direction of the stream at infinity. 



The Lengths of the Planes. 



The lengths of the planes are given 



by I 



dz 
dt 



dt and 



i 



dz 



dt 



dt. Considering 



/*i 



the latter, it can be written in the form 



2ir 



B 



dz 



dJJ, or, since = along the rear 



'irr 



plane and z real, in the form dz, z^ and za„ denoting the values of z when 



ZB 



U = B and Jtt respectively. 

Now f = dz/dw, and to = (sin^ U — sin^ P)^ ; 
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therefore 

L2 (the length of rear plane) = 4 ^ sin U cos U sin (U -j- P) sin (U — P) ^U. 



'•Itt 



B 



Substituting for ^ from (4) log ^ = O, it is immediately reduced to the 
form 



L3 = 2 

Similarly Li = 2 
for the front plane. 



sm 

B 

rB 



in2(lT + P)sin2U i^^^Lg^^JT'^'' ^U. 

Lsm(U — B)J 



.sin2(U + P)sin2U /^!^^^"^ST'^''^U, 
^ ^ \sm(B-U)J 



(8) 



(9) 



The Pressures on the Flames. 



The pressures are given by 
•^ dz 



iJ<'-«"> 



dt 



dt and ^ (1—$^) 



dz 
~dt 



dt. p = density of fluid. 



We first of all proceed to find 



>2 



dz 
dt 



dt for the rear plane. This and the 



corresponding integrals will be called the Pressure Difference Integrals. 
^ =r 0, as before, along the rear plane, hence we write 

fdwY 7 ^dw 

Proceeding as in the case of the lengths, we obtain 

= 2 f ^in2 (U - P) sin 2 U I ^!"^,^~^J T "^"d^, 
J„ ^ ^ Lsin(U + B)J 



P2 = 


q^ dz = 


• 


2b J 



and, similarly, 



Pi= 2 



Iin2 (U - P) sin 2 U I ^-"S , ?t\ \ ^'^"'^'V- 



(10) 



(11) 



§3.-4 Further Amplication to the Case of a Flane with more than One Bend. 

Suppose now that we have a plane bent at several points A2, Ag..., A,...., 
An~\. Let the angles the parts A1A2, A2Ay, Aj^A^+i..., A^-iA^, make with the 
direction of motion of the fluid at infinity be ^2, «n «r+i; •.., ^n respectively. 

Let the point P where the stream divides be situated between A^ and 
Aj.+i, and the points of maximum velocities be E2..., Er, ^\, E^+i..., 
respectively. 

AiE and A^B are the free stream-lines along which ^ = 1 as before 

The z, X2, w, and t diagrams are as shown in figs. 6, 7, 8, and 9. 
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Ai 
Fig. 8. 






E 


Ai As ^i^ A^ 


P 

1 


Ar+i 


A„_i 


Aji 



E at —00 



Fig. 9» 



E 



D at 00 



r+l 



We now apply the Schwarz-Christoffel formtila to transform fig. 7 to fig. 9 ; 
in this case it takes the form 

dXl _ const, {t — p2) (t — pr) (t — Pr) (^ — pr+l) 

dt {t-—a2){t--ar){t'-f){t-—ar^\){t'—an-i)^/{t'-'ai){t--'aj^ ' 
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p2, pr, p^r, pr+1 coriespond to the values of t at E2, Rr? RV> Rr+i, and the 
a's to the values t corresponding to the A's. 

For more bends in the plane additional factors would be inserted in the 
numerator and denominator of this expression. 

How, putting as— 5^2 = 7r/2r2, ar+i — ccr = ir/2vr, etc., we have, on 
splitting up the right-hand side of this equation by means of partial fractions 
and determining the constants as in the previous case, 

dO^ _ x/{ap'—ar){ctp'—cti) , 1 \/{ci2 — ctn) {0-2 — ^1) ■ 

. 1 v^((^y — an) (% -— ^^1) , _.!__ \/(^r4-l — %)(<^yi-l — %) ,_ 

2vr ' (t--ar)^/{t—ar)(t—ai) 2vr+i ' it--ar+i)y^{t—an){t—ai) 

I 1 vC^yt-l — <^'»)(%-l — %) (12) 

2vn~i ' {t — an-i)\/{t—a^){t — aiy 
Integrating as before, and remembering that ^ = e^, we have 

f = J'\ /{^n—ctp){t—al)-\-^JXal — ap){t — an) 
L \/{an—cii){t—ap) 

To find the angle of attack on the part An+iA« (fig. 6), put ^ = oc in the 
equation for XI (not given here, but obtained from (12) by integration). 

O __ 2 loo- \/(^^~^^y)"^ \^((^l^Ct)p , V A JQo- yy^n — ^r) + \/{ch — Ctr) 

Kow, % > %-i > . . . > a^+i > ^r > . . . > ch, hence the imaginary part of 
this expression 

or using the same artifice as before (3, § 2) 

= i(2P+ S iu. , 
where TJy is the value of U corresponding to a^. 

r — n—l 1 

Therefore = 2P+ 2 - U^ (compare 6, |2). 

r = 2 I'r 

6 is known from the conditions of motion about any particular plane, 
U2, ..., Un-i are also known (depending on the points at which the section of 
the plane is bent) ; hence P can be determined. 
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§ 4. Reduction of the Integrals, 

When l/2i/ is an exact fraction Ijn say, the integrals are reducible to finite 
form. ^, 

Taking, for example, the integral 

psin2(U + P) sin 2U i^^^2-Sr^V^''^U, 
Jb lsm(U — B)J 

, ... sin(U + B) 

and wntmg . )^^ — =r{ = z^, 

sin(U— B) 

it will be seen, by expanding the first two factors of the integrand in powers 
of sin U and cos U, that the integral can be expressed as the sum of integrals 
of the form 

sin^ tJ cos^ TJzHJ] , where ^ + (;^ = 4. 



1 



If we now change the variable from U to z, we have 

tan \J __z^-\-l ,jj _ —nz^~'^ sin 2^djZ 

tanB ~"^^-l' ^ "";^2«_2;^n(3og2B + l' 

TT _ (;2;'* — l)cosB . • tt _ (^''* 4-1) sin B 

cos U — — , . o — 7z~^ ttj anQ Sm U — — r- — -z — a:"-^ tt* 

y^(;^2n_2;S« cos 2B + 1) ^y{Z^''-2z'' COS 2B+ 1) 

The integrals are all thus reduced to those of rational fractions having as 
denominator (2;^^ — 22;'^cos2B-f 1)^ and as numerator integral powers of z^ 
and these can easily be integrated by means of partial fractions. While the 
functions involved have been shown to be " integrable," the results are 
obviously far too cumbersome to be suitable for the purpose of calculation, 
and we therefore have recourse to approximate methods of integration. 

§ 5. ApjfToximate Uvahiation of Integrals. 

'- Simpson's Eules " will be used to evaluate the integrals where they are 
applicable, but it is evident that in the neighbourhood of U = B, they fail. 
The following method will consequently be used near U = B. 

Consider sin^ ( U ■— P) sin 2 U < ^]^ ;^^'^ ^{ r dXJ, 

J Lsm(U + B)J 

Put U— B = a^, and the integrand becomes 



sm"* X 



sm^(a^ + B--P)sm2(a? + B) -r./ x 
^i : ^ ^ — 1 — i = jj (x), say. 



When x—*0 (i.e. U —* B), sin*^ x —* x'"\ and we may write F (x) = x^^ (x). 

(j) (x) is now expanded in powers of x, being finite and continuous within 
the limits, and F (x) = (A^*^* + B^*^+i + Gx'^'-^^ +...) = 2/- 
Let yi be the value of ^ [ = F (x)] when x = h, 

« 
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rrh 

Suppose yclx (to be required) = (P13/1 + . . . + ^rVr) ^^' 

Jo 

= (Pi + 2"T2 + . . . + r^^Pr) AA^'^+1 + (Pi + 2"^+iI^2 + . . . + r^-^^^r) B/i^^^ ^^^^ 
But ^ydm = A (r/^^"^ V(m + 1) + B (rhy'^^/{m + 2)4-.... 



Therefore Pi + 2^T2 + S^Tg + . . . + rH\ = r«*+7(f/i + 1), 

Pj4.2«^+iP2+3"^+iP3+...+^^'''"^^Pr = r^+V(^ + 2X and so on. 
which, when solved, give the values of P. 

In the cases that will be given the lengths and pressures have been 

calculated for values of m = 1/12, 1/18, 1/36, 1/180, and 0, corresponding to 

angles 7r/2p = 15°, 10^ 5^, 1°, and 0° between the planes, and r has been 

taken = 2 giving then 

Pi = 2«^+7(m4-l)(m + 2), 

P2 = 2 m/(m + 1) (m -f 2). 

The formula here derived, then, will be used in the immediate vicinity of 
B, from B to Bi say, on the front plane ; 0<Bi<B ; and from B to B2 on the 
rear plane, B<B2<-|7r. One of Simpson's Eules will be used for the 
evaluation of the integrals between Bi and 0, and B2 and -^tt, the intervals 
being 5° 

To test the accuracy of this procedure, the integral .. g" .g. has been 

Jo\/\^ — -^ } 

■evaluated and found to be 1*55, as compared with the true value |-7r = 1*57. 
The error is, therefore, about 1*5 per cent, with these intervals. 

In applying our rules to the integrals, there are two ways of proceeding. 

(a) Certain fixed values can be given to P and B, P having, of course, to 
satisfy the condition assigned to it in (7), § 2, and the total lift and drift 
with the angles of attack determined by these values of P and B found. 

(b) A fixed value can be assigned to B w^hich will determine the ratio of 
the length of the front to the rear plane, and we have three integrals 



sin2U.Z«^^^U, 



sin^2U.Z^^rro, 



sin4U.Z-t^lJ, 



to evaluate, with the coefficients |, | sin2P, and — |cos2P respectively, for 
half the lengths of the planes. For sin^ (U + P) sin 2 U can be written in the 

form 

J sin 2 U + i sin 2 P . sin^ 2 U — | cos 2 P . sin 4 U. 

[In the pressure difference integrals Z is different, and the sign of P is 
changed.] 

This method has the advantage that after the integrals have been evaluated 
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we can find the thrusts for any angle of attack, which will determine for u.s 
the value of P. One or two cases will he dealt with according to the first 
method, but the second way will in general, in what follows, be used. 

Now, after the lengths of the planes for various values of B have been 
calculated, they will be reduced in such a way as to make the sum of the 
semi-lengths equal to 10 e, where e is the circular measure of 5° 

The corresponding lifts and drifts for various angles of attack, and values 
of 7u, will accordingly be calculated. Further the thrust on the chord AG' 
(fig. 1, I 2) will be calcuilated, and results compared. 

Finally the lifts and drifts on planes, for different values of B, and the least 
integral value of P satisfying the condition already cited, will be calculated. 
This will be done for one value only of m, this being sufficient to give a general 
conception of the relative advantages of the planes considered. 

Calculations have been made for the following values of B ; 35°, 40°, 45°, 
60°, and 75° and for the values of m already given. For any other values at 
B and m the corresponding results can be found by interpolation. 

The thrust on the chord is easily obtained as follows : — 

A 



N 




M 



Fm. 10. 



BN and CM are parallel to the direction of motion, AHM and BP perpen-^ 
dicular from A and B (fig. 10). 

The angles of attack ABN and BCM are a and ^8 respectively. 

Let AB = li and BC = k (which are supposed to be already calculated). 

AM = AH + NM = 4sina + 4sin^, 

CM = CP + PM = Zi cos a + /s cos /3. 

Therefore the angle of attack (f = ACM) on the chord is 

tan"i ^1 -"'"^ ^ + ^2 sin ^ 
li cos « + 4 cos ^ ' 

Knowing this, P can be found, and then by means of the method (6) above 
with m = the thrust can be found. 
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§ 6. N'lmnerical Results. 
The following Tables are given, to illustrate the way in which the general 



results were obtained. 

In the first case (Table I) — 

B = 45°. P = 5°. 
Length of front plane, 2-78 x 2 c. 
„ rear „ , 7*22 x 2 c. 



9 = 

a. = 

Q = 



circular measure of 5°. 

density of fluid. 

angle of attack on front plane. 

,, jj rear „ 

velocity of fluid at infinity. 



Table 1. 



m. 



/3 — a. 



Front. 



Rear, 



a. 



L/pcq' 



D/pcQl 



X 

1 '2 
_1_ 

18 

JL. 

30 

18 





o 


o 






15 


2i 


1-638 


0-068 


10 


5 


1-361 


0-118 


5 


7k 


1-351 


0-177 


1 


9i 


1-216 


0-202 





10 


1-188 


0-208 



/3. 


L/pcQ-l 


B/pcQ2. 


o 

m 


2-230 


0-702 


15 


1-912 


1 0-512 


12^ 


1-568 


1 0-346 


lOi 


1-268 


1 0-233 


10 


1-188 


! 0-208 



Sum. 


L/pcQl ' B/pcCll 


3 -868 
3-273 
2-919 

2-484 
2-376 


0-770 
0-630 
-523 
0-435 
0-416 ! 



Chord. 



L/pcQ,". 



B/pcQl 



2-976 
2-778 
2-678 
2-410 
2-380 



0-698 
-554 
0-530 
0-432 
-414 



The ratio of the lengths of the two planes has been found in all cases to be 
approximately independent of m and P. It does not deviate more than 
1*5 per cent, from the mean value. 

The next Table (II) gives the total lift and drift on a plane in which B = 75° 
for the cases (1) P = 2°, and (2) when the angle of attack on the front plane 
is zero. The lift and drift on the chord are also given. 

The lengths of the front and rear planes in this case are 8*77 x 2g and 
1-23 X 2e respectively. 

Table II. 




Bent plane P = 2°. 



L/p.'Q2. 



D/pcQl 



Bent plane. Front attack 



LjpcQ;. 



'D/pcQ:' 




2-434 
1-918 
1-479 

1 -228 
1*020 



0-150 
0-113 
0-092 
0-080 
0-070 



2-111 


0-084 


0-953 


0-057 


1 -388 


-036 


0-946 


-058 


0-687 


0-010 


0-958 


0-062 


0-189 


Negligible 


1-006 


0-069 








1-020 


0-070 
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Thirdly, the lift and drift have been calculated for the least integral value of 
P satisfying the conditions of the problem, and the results given in Table III 
m = 1/36, /3— -a = 5°. 

Table III, 
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^ sum of 


B. 


P. 
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40 
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2-01 
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0-923 


1-833 
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2-756 


0-412 


45 


4 


2-78 


7-22 
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1-394 


0-111 


0-286 


2-551 


0-367 


60 


3 


5-84 


4-26 


1-655 


-487 


0-125 


-088 


2 -132 


0-213 


75 


2 


8-77 


1-23 


1-340 


0*139 


0-073 


0-019 


1-479 


092 



Finally in Table IV are given the total lift and drift, and the ratio of lift 
to drift, for a plane whose angle of attack on its front part is 5°, for values of 
B = 35°, 45°, 60°, and 75° and m = 1/12, 1/18, 1/36, 1/180, and 0. 

Table IV. 



B. 


35°. 


46^ 


60°. 


75°. 


fi — a. 
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4 -182 
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1-544 
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0-154 
-114 


3-41 

4-61 

6-55 
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1-016 
0-628 
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5-36 

7-05 

10-36 
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3-970 
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2-315 
1*551 
1 -285 


0-690 
-462 
0-269 
0-148 
0-113 


5-76 

7-12 

8-94 

10-50 

11-45 


3-110 
2-509 
1-891 
1-462 
l-2b2 


0-366 
0-263 
0-178 
0-130 
0-112 


8-53 

9-54 

10-62 

11-25 

11-45 



The results contained in this last Table are also given graphically as curves 
where the lift and drift are plotted against ^ — oc, fig. 11. The curves 1, 2, 3, 
and 4 represent the lift for B = 35°, 45^ 60°, and 75° respectively, while 
(1), (2), (3), (4) represent the drift. 

A curve is also given, in which the length of the front plane is plotted 
against values of B (fig. 12). 

In conclusion we now proceed to examine the above Tables. It is evident 
that in all cases both the lift and drift are largest when the angle between 
the planes is largest (for the values of the angle under consideration). This 
evidently points to the fact that the ratio of lift to the area of the plane 
is greater, the greater the bend. But; the drift is also increased, which is 
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not desirable in the case of an aeroplane. Again the ratio of lift to drift is 
greatest when the angle between the planes is least, but the lift is small. 

These two conclusions agree with the experimental results that aerofoils 
with a large camber give a better ratio of lift to surface, whereas nearly flat 
aerofoils give a better ratio of lift to drift. 

The figures given in Table I show, however, that the ratio of lift to drift, 
for a large value of m even, is greater than the corresponding ratio for the 
chord, and the lift on a bent plane is much greater than on its chord. Hence 
it appears that a bent plane is decidedly more efficient than its chord. 

Further, from the last Table it is seen that as the length of the front plane 
increases, the total lift decreases, whereas the ratio of lift to drift increases. 
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Hence a plane bent near the leading edge gives a better ratio of lift to 
surface, while on the other hand a plane bent near the rear edge gives a 
better ratio of lift to drift. Combining these last conclusions with the 
corresponding ones above, a plane with a large bend close to the front edge 
gives the best ratio of lift to surface, while a fairly flat plane bent near the 
rear edge gives the best ratio of lift to drift. 

All these conclusions, then, tend to show that a distinct advantage is gained 
by making use of a cambered plane. 

The above methods can be applied to find the thrusts on any bent plane if 
the ratio of the lengths of its front and rear parts be given. Fig. 12 gives 
immediately the corresponding value of B. Having been given the angle of 
attack as well, and the angle at which the plane is bent, P and m are easily 
found. The case of a plane bent twice can also be considered (as shown 
in § 3) but the integrals become complicated. 

The following application of fig. 12 may, in certain cases, be of use. A glance 
at fig. 12 shows that between certain limits (viz., between 37° and 75*^ approxi- 
mately) the leugth I of the front plane is a linear function of B = aB + Z>, 
say. If, then, the values of B determining the length of the middle part of 
the plane lie between these limits, we may proceed thus. 

Equation 8, § 2, gives in this case 



z 



2 sm^(U+P)sm2U^ . ).^ — ^^ i . ;^^ — -{V dJJ 
Jb, Lsm(ir~Bi)J Lsin(B2-U)J 



where B2 > U > Bi, B^^ and B2 being the two values of B, 

and vi and V2 are the angles of bending at Bi and B2 , 



IBs 

aU + h 



therefore 

a - 



2sinHU + P)sin2u|?|iL(J^^ 

^ ^ ^ Lsiii(U~Bi)J Un(B2~-U)J ' 



^^ r sin(U + Bin 1/2-^ rsiir(B24jJ)"l ^/^.^ ^ g 

Lsin(U-Bi)J lsin(B2-U)J 2sin2(U + P)sin2U^ 

which, when substituted in the formula corresponding to 10, § 2, gives the 

pressure difference 

4 f B2 

- sin2 (U + P) sin2 (U - P) sin^ 2 \Jd\J, 

It should be observed that this is not independent of vi and V2, these being 
included in P. This integral reduces to 

- \ sin2 2 U sin* U c?U - ^ sin^ P f sin^ 2 U sin^ I J ^U -f - sin* P [ sin^ 2 U ^U, 
a} cc J <^ J 

the three integrals being easily reducible.* 

* Williamson, * Int. Calc.,' Ch. III. 
VOL. XCI. — A. 2 G 
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This metbod has been tested and found to give results which agree closely 
with those obtained by the more laborious method applied to the case of the 
plane with one bend. 



The Difference between the Magnetic Diurnal Variations on 
Ordinary and Quiet Days at Kew Observatory, 

By C. Chree, Sc.D., LL.D., E.K.S. 

(Eeceived March 24, 1915.) 

§ 1. While investigating the nature of the diurnal inequality of the magnetic 
elements at Kew Observatory, I have discovered a somewhat remarkable 
relationship between the inequalities derived respectively from quiet days 
and from ordinary days {i.e. all days except those of large disturbance). 
The period whose records were employed consisted of the 11 years 1890 to 
1900. During that period five days were selected for each month by the 
Astronomer Eoyal as typical of quiet conditions. An analysis of the diurnal 
inequalities at Kew on these quiet days was given in an earlier paper.* Use 
is also made of the diurnal inequalitiesf derived from 209 disturbed days, 
being the days of principal disturbance from 1890 to 1900. The new 
material employed consists of diurnal inequalities based on all the days of 
the 11 years, excluding the 209 highly disturbed days but including the 
quiet days. The cost of measuring the curves was defrayed by a grant from 
the Government Grant Committee. 

The diurnal variation in the horizontal components of magnetic force is 
often advantageously studied by considering the vector which represents the 
force to whose action the departure from the mean value for the day may be 
ascribed. 

In the figure NS and WE are perpendicular lines, respectively in and 
perpendicular to the geographical meridian. If CM and PM, the 
co-ordinates of the point P, are respectively proportional to the departures of 
the east and north components of force at any given hour from their mean 
values for the day, OP represents in magnitude and direction the force to 
which the diurnal variation in the horizontal plane may be ascribed. The 
successive positions occupied by P throughout the 24 hours constitute what 
is known as the vector diagram. 

* ' Phil. Trans./ A, vol. 202, p. 335. 
t * Phil. Trans., A, vol. 210, p. 210. 



